Semiclassical dynamics of atomic Bose-Einstein condensates 
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An atomic Bose-Einstein condensate (BEC) is often described as a macroscopic object which can 
be approximated by a coherent state. This, on the surface, would appear to indicate that its behavior 
should be close to being classical. In this paper, we clarify the extent of how "classical" a BEC is 
by exploring the semiclassical equations for BECs under the mean field Gaussian approximation. 
Such equations describe the dynamics of a condensate in the classical limit in terms of the variables 
{x) and (p) as well as their respective variances. We compare the semiclassical solution with the 
full quantum solution based on the Gross-Pitaevskii Equation (CPE) and find that the interatomic 
interactions which generate nonlinearity make the system less "classical." On the other hand, many 
qualitative features are captured by the semiclassical equations, and the equations to be solved are 
far less computationally intensive than solving the CPE which make them ideal for providing quick 
diagnostics, and for obtaining new intuitive insight. 



PACS numbers: 03.75.-b,03.75.Kk,03.65.Sq 
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An atomic Bose-Einstein condensate (BEC) is a 
macroscopic quantum object composed of typically mil- 
lions of bosonic atoms. By now it is well-established that 
a large class of phenomena involving atomic BECs can 
be described accurately by the Gross-Pitaevskii Equa- 
tion (GPE) for the condensate mean fieldd. The inco- 
herent component outside the mean field can, in most 
situations, be ignored as they become significant only at 
higher temperatures closer to the critical temperature of 
BEC phase transition^. The GPE, which is a nonlinear 
Schrodinger Equation for the condensate wave function 
ip in the presence of an external potential V ex t(x), is 
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where Uq = A-kH aNy/m with a being the s-wave scat- 
tering length; N is the total number of atoms in the 
condensate. 

To derive the GPE, the quantum state of a BEC is 
assumed to well-approximated by the coherent state 
the best justification for this choice of quantum state is 
of course, the success of the GPE itself. There are, in 
fact, a number of unique characteristics associated with 
the coherent state, and the one that has gained most at- 
tention in the context of BEC studies is that of the phase 
coherence - an indication that the many-body quantum 
system is acting effectively as a single quantum object. 
Another important characteristic of the coherent state is 
that it is the quantum state closest to the corresponding 
classical system. A coherent state \a) can be shown to 
be identical to minimum uncertainty Gaussians where: 
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and it was Schrodinger himself who first showed that 
Gaussian wave packets form coherent states of a har- 
monic oscillator that display classical motion. Further 
studies in various contexts [I, H, 0] have confirmed that 
Gaussian dynamics is indeed classical at all levels H Q. 



We note that the "Gaussian dynamics" in this context 
includes a "squeezed coherent state" i.e. Gaussian of 
variable width, and the state is assumed to remain a 
Gaussian throughout the evolution. 

Despite the coherent state/Gaussian ansatz, the solu- 
tion to the GPE does not necessarily give classical dy- 
namics. This apparent "contradiction" arises simply be- 
cause the coherent state assumption used to derive the 
GPE is made at the level of second quantization - the 
ansatz turns the many-body quantum field operator into 
a classical field of condensate wave function. The result- 
ing wave function in real space can, in principle, take 
any form depending on the boundary and external con- 
ditions. In order for the classical dynamics to emerge 
from the condensate wave function, a further "coherent 
state" assumption within the mean-field level of theory - 
that the condensate wave function takes on a Gaussian 
profile - is needed. The fact that a BEC is not strictly 
a "coherent state" within the mean-field regime is evi- 
dent from examining the Q-function [|| for the conden- 
sate at the level of GPE: a Gaussian wave packet gives 
a circular Q-function expected of a coherent state while 
the Thomas-Fermi approximation for which most exper- 
iments are valid, gives an elliptical Q-function indicative 
of its nonclassical nature @. 

The Gaussian profile approximation can, in fact, be a 
good approximation for the stationary condensates with 
a small number of atoms or with an attractive inter- 
atomic interaction. However, the inherent nonlinearity 
of the system can destroy any Gaussian profile approxi- 
mation very quickly during the dynamical evolution. We 
note that a Gaussian description was used previously in 
the context of studying surface collective mode of a BEC 



with attractive interactions 1C | , where the results based 
on variational methodsfHI, 12 were used. In the early 



works based on variational principle [Til [l2| . it was shown 
that many of the features of a BEC can be captured using 
the Gaussian ansatz. However, the semiclassical aspect 
and the validity of the Gaussian dynamics compared to 
the full quantum dynamics was not fully explored. Also 
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the results were restricted to the harmonic trapping po- 
tential, which doesn't demonstrate a more general picture 
of deviation from the GPE. 

In this paper, we ask the question "Just how classical 
is a BEC and under what conditions?" We explore the 
dynamics of the atomic BECs in the semiclassical limit 
and compare with the full quantum solution described by 
the GPE. The semiclassical limit is obtained by assum- 
ing that a BEC maintains a Gaussian profile (of variable 
mean and variance) throughout the evolution and de- 
riving the corresponding effective Gaussian semiclassical 
equations. Strictly speaking, only the ground state for 
a harmonic trap with Uo = in the GPE of Eq. ([I]) is 
precisely Gaussian; for Uo ^ 0, a Gaussian profile pro- 
vides an idealized classical limit that acts as a "baseline" 
for comparison with the full quantum result. The mo- 
tivation for this study comes not just from addressing 
fundamental questions on the quantum-classical nature 
of a BEC but also in the context of advancing quantum 
engineering and control methods [Hj via simpler equa- 
tions that relate relevant observables. The semiclassical 
equations can provide us a means to understand complex 
quantum behavior in terms of the more familiar and in- 
tuitive classical picture, providing a new insight into how 
to manipulate atomic BECs. 

The usual derivation of Dirac's time-dependent vari- 
ational principle states that an action of the form 
r = / dt(V,t\ih(d/dt) - H\V,t) 5T = results in the 
Schrodinger Equation. The solution is obtained by re- 
stricting \^>,t) to a subspace of the full Hilbert space. 
For deriving the semiclassical Gaussian dynamics, the 
subspace is assumed to be that of a "squeezed state" 
wave function of the formQ: 



ip(x,t) = 



exp { -^(a; - (x)f + ip(x - (x)) 
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where, with the definition AA = A — (A), a = 1 
i(AxAp + ApAx) and p is defined such that 



(Ax 



2m\ 



p 2m (2m\) 



(3) 



with m = 1 corresponding to the position variance, 
(Ax 2 ) = p 2 . The use of a squeezed state ansatz is rem- 
iniscent of the time-dependent Hartree-Fock-Bogoliubov 
(TDHFB) theory [14J . The major difference is that with 
TDHFB the squeezed state ansatz is applied for both 
the condensate and non-condensate atoms, as it is a fi- 
nite temperature theory. Further, the mean field solu- 
tion used is the precise ground state for the system, not 
a Gaussian approximation i.e. it has nothing to do with 
the classical dynamics. TDHFB is also a very difficult 
problem to solve. Here, we are interested in the semi- 
classical dynamics at the mean field level, i.e. where all 
the atoms are assumed to be in the condensate. 

With such an ansatz, the dynamics of the system char- 
acterized by a Hamiltonian H = p 2 /2 + V(x) (here we 
set the mass to unity) can be represented as an extended 



classical gradient system for the average and fluctuation 
variables. The dynamical equations for the coordinate 
variable (x) and its standard deviation, p areQ 
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where V^((x))= £^V{x)\ x=(x) . 

By integrating the above equations once with respect 
to time, the dynamics of the respective conjugate vari- 
ables (p) and n are obtained: 



where 
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From the uncertainty relation 
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the momentum variance can be written in terms of p and 
n as 
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implying that whenever n = 0, we have the minimum un- 
certainty state. We note that Eqs. (|H - E|) are the most 
general equations of motion for (x), (p), p and n for a 
squeezed coherent state in an arbitrary potential V(x). 
We use this set of self-contained equations to model our 
case. It is noted that these equations may also be derived 
alternatively using an operator expansion method[l5[. A 
point of interest is how to include the interatomic inter- 
action in the equations of motion. 

As shown in Eq. {T]), the GPE for the condensate 
wave function is a nonlinear Schrodinger equation with 
an additional mean field potential term Uq\i(j(x, t)\ 2 . We 
therefore model nonlincarity by extending our potential 
to introduce the "internal" mean field potential expe- 
rienced by an atom in the condensate i.e. V(x, t) = 
V ex t(x,t) + Vi n t(x,t) where V ex t(x,t) is the usual exter- 
nally imposed potential (such as the harmonic trap) and 
Vi n t (2, t) is the internal "mean field" potential due to the 
interatomic interactions experienced by the condensate 
atoms: 



V int (x,t) = U \^(x,t)\ 2 = 
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where we have assumed that the condensate wave func- 
tion is given by the squeezed state Eq. ^j. Since the 



3 



semiclassical equations are essentially equations govern- 
ing (x) and p and Vi nt (x, t) is itself a function of the time- 
dependent variables (x) and p, it should provides the nec- 
essary nonlinear effect consistent within this framework. 
We shall check this by compa ring with the equations de- 
rived independently in Ref. pjj where the Lagrangian 
itself contained the nonlinear interaction term from the 
outset. 

For the purposes of deriving the semiclassical equa- 
tions for atomic BECs, we shall consider a most general 
(practical) form for the external potential: 



V ext {x,t) 
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a[cos(fcx - nt) -l]+0x (11) 



where to provides the frequency of the harmonic trap, a is 
the amplitude of an optical lattice, [3 is the gradient of the 
tilt of the optical lattice, and O allows the possibility of a 
travelling lattice at speed fi/fc. By setting a = f3 = CI = 
and uj ^ one is considering the case of a harmonically 
trapped condensate, while setting a ^ and 10 = would 
imply the presence of an optical lattice only. It is noted 
that all of these parameters a, /?, k, fl and to can, in 
general, be time-dependent quantities. The form of the 
optical lattice potential was chosen so that at time t = 
it approximates the standard harmonic potential for a 
condensate initially at mean position (x) = 0. 

It was found that the infinite series of Eqs. ([4]) and (O 
involving the V ex t.( x ^) part can be obtained relatively 
straightforwardly, while that involving Vi n t{x,t) was less 
straightforward, and required the use of the definition of 
the Hermite polynomials H n {x) 1 



H n (x) = (-l)»e* — e- , 
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the identity 

H 2n +i(x) = (-l)"2 2 "+ 1 n!x4 1 / 2 )(x 2 ), (13) 



where L„ (x) are the generalized Laguerre polynomials, 
and the generating function for the Laguerre polynomials 
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After some work, the semiclassical equations that de- 
scribe a BEC in the general external potential and in 
the presence of a mean field is finally obtained as: 
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These equations in (x) and p can be very easily solved 
numerically compared to the GPE which requires a dense 



grid in both position and time for accuracy. In particular, 
the equations are essentially identical to those obtained 
in Ref. This proves the validity of the assumption 

that at least for the Gaussian anstz, the effect of non- 
linearity can be modeled exactly by assuming that the 
interatomic interaction effectively results in a modified 
potential. This insight can help in analyzing equations 
in other contexts. 

The variables in the semiclassical equations yield an 
extended potential system where the fluctuation and av- 
erage variables are treated on equal footing. The ex- 
tended Hamiltonian H ext . = (H) in our case is given by 
H ext . = H x + H p + H xp where 



M! +/ 3(x) + iu, 2 (x> 2 (17) 

-acos(fc(x) ~nt)e- p2k2/2 (18) 

n 2 i 2 2 h 2 u 

-Y + r p + w 2+ ^b- (19) 



The form of the extended Hamiltonian indicates effec- 
tively an infinite barrier potential H ex t. — > oo as p — * 
i.e. the "quantum fluctuations" can never go to zero ex- 
cept in the limit h — * and Uq — > 0. The last two terms 
of Eqs. (|16p and (|T9|) indicate that the effect of having a 
mean field interatomic interaction can be thought of as 
effectively modifying the Plank constant h — > h e f /. such 
that: 
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This formally supports the observation that nonlinear in- 
teratomic interactions (having Uq > 0) push the system 
deeper into the "quantum" regime. 

Several interesting observations can be made simply 
by examining Eqs. (fl5|) and (fT6|) a bit more closely. For 
instance, it is clear how to directly set the accelerations, 
d 2 {x)/dt 2 and d 2 p/dt 2 by choosing various initial param- 
eters a, k, f2, ui, (3 and Uq. This can be useful in designing 
quantum control strategies to manipulate a wave packet's 
acceleration. Such direct and intuitive information is not 
readily available from the GPE. In particular, the entire 
right hand side of Eq. (fTT)|) can be made to vanish i.e. 
one can figure out how to stop the variance of a wave 
packet from accelerating. 

This is significant since, if one starts off with a station- 
ary state i.e. with a zero initial rate of change dp/dt = 0, 
the width of the wave packet throughout the evolution is 
"frozen" at the initial value under such conditions. For a 
Gaussian this is tantamount to having its overall shape 
unaltered throughout its evolution. Although such a be- 
havior is reminiscent of a soliton solution [17] . they are 
better classified as solitary waves [18| as it has less to do 
with the nonlinearity of the system and more to do with 
being a stationary state of the extended Hamiltonian. 
Another way to look at this is to consider the gradient 
of the extended potential with respect to the position 
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FIG. 1: Solid line represents the GPE solution and dashed 
line represents the Gaussian semiclassical prediction, (a) Po- 
sition probability density |f/>(a;)| 2 for the "stationary state" 
corresponding to Uo = 5. (b) the width of the stationary 
state wave packet represented by the position variance for 
various values of Uo- (c) evolution of the position variance of 
the Uo = 5 wave packet after being released from the trap, 
(d) evolution of the momentum variance of the Uo = 5 wave 
packet after being released from the trap. Time is in units of 
harmonic oscillator period throughout this paper. 
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p continues on with its initial value without acceleration. 

For the case of a harmonically trapped condensate 
(a = f3 = 0, u) > 0) the right hand side of Eq. (fTB)) van- 
ishes when the initial width of the wave packet is such 
that 



P = 



J eff. 
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i.e. it gives the expected width of the stationary state 
solution for a given value of Uq. In Fig. Q] the predic- 
tion of the semiclassical equation for the stationary state 
solution is compared with those obtained u sing the full 
quantum ground state solution of the GPE[lg. Figure 
[IJa)compares the shape of the Gaussian solitary wave so- 
lution for Uo = 5 with the ground state of GPE for the 
same value of nonlinearity, while Fig. [IJb) compares the 
width of the stationary state wave function as a func- 
tion of Uq. We have checked using the numerical ground 
state solution to GPE with a harmonic trap[l6| that for 
Uq up to around 9, the profile of the condensate can be 
well-approximated by a Gaussian before going over to the 
inverted parabolic form of the Thomas-Fermi regime. 

A common experimentally relevant situation is when 
the harmonic trap is turned off in the absence of an op- 
tical lattice, (q = (3 = 0, uj = 0, Uo > 0). It is easy 



to see from the right hand side of Eq. (| 16|) being pos- 
itive that the position variance p will start to increase 
without bound when the trap is suddenly turned off. In 
addition, as p — > oo, d 2 p/dt 2 = dIL/dt — > due to the 
last two terms of Eq. (p~6|) i.e. II becomes a constant, 
and consequently 
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(23) 



i.e. it is expected that the momentum variance reaches 
a constant value after a while. Although this fact about 
the momentum variance is not at all obvious from the 
GPE or even from the theory of Fourier transforms - 
one would suspect the momentum variance tends towards 
zero as the position variance increases without bound 
Fig. [ljc) and (d) indicate that the quantum simulation 
using GPE shows a good qualitative agreement with the 
semiclassical prediction. 

What is also notable is that, if one wants to hold the 
wave packet together in the absence of any confining po- 
tential (a — (3 = 0, ui = 0), one needs Uq = —H 2 ^/tt < 
to get the right hand side of Eq. (jTHJ) to vanish i.e. an at- 
tractive interaction, as expected. This also presents the 
limit as to how far negative Uq can be for a harmonically 
trapped BEC with an attractive interaction. This is be- 
cause, with a harmonic potential present (a = (3 = 0, 
uj > 0), Uo = —h 2 y/TT in Eq. (p~6|) eliminates the infinite 
barrier potential so that p can take on zero value after 
some time. This then makes (Ap 2 ) = |- +LT 2 to diverge, 
making the system unstable - this limit on Uo shows that, 
according to a semiclassical argument, a given harmoni- 
cally trapped BEC with a negative scattering length can 
support up to 



N = 



1 



0.141 



4^|c 



(24) 



atoms where the interatomic scattering length a is in 
harmonic oscillator lengths yKfmw. Ignoring additional 
considerations such as the dimensionality of the conden- 
sate this is of comparable order of magnitude to the pre- 
diction made using the full quantum theory [l!| for the 
upper limit in atom number, N — 0.573 /j aj . Again the 
semiclassical equations can be a useful tool in finding 
order-of-magnitude predictions based on intuitive argu- 
ment. 

In a previous work[20j, the case where the interaction 
strength Uo of a harmonically trapped atomic BEC is 
suddenly tuned to zero was studied using the full quan- 
tum theory. In that work it was noted that a BEC starts 
to undergo oscillatory motion which correspond to the 
rotation of the error contour. We investigate here a simi- 
lar effect using the semiclassical equations: starting with 
the stationary state with nonlinearity Uo > 0, the non- 
linearity is suddenly switched to U new at time t — and 
the subsequent evolution of the position and momentum 
variances is observed. We consider four possible cases: 
(a) U new = 0, (b) < U new < U , (c) U new > U and 
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FIG. 2: Time evolution of (Ax 2 ) (Solid line for the GPE sim- 
ulation and dash-dot lines for the semiclassical prediction), 
and (Ap 2 ) (dashed line for the GPE simulation and dotted 
line for the semiclassical prediction) expected when a BEC 
with the initial nonlinearity of Uq — 5 is tuned to new non- 
linearity U new at time t = 0. (a) U ne w = 0, (b) U new = 2, (c) 
U n ew = 7, and (d) U ne w = —0.5 
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FIG. 3: Time evolution of (Ax 2 ) (Solid line for the GPE 
simulation and dash-dot lines for the semiclassical prediction) , 
and (Ap 2 ) (Dashed line for the GPE simulation and dotted 
line for the semiclassical prediction) expected when a BEC 
with the initial nonlinearity Uq = 5 and position variance p 2 is 
placed in an optical lattice of different wavelengths, Xol- (a) 
Xol = p, (b) Xol = 2p, (c) Xql = 10p, and (d) Xql = 20p. 



(d) —h 2 y/ir < U new < and compare the quantum cal- 
culations with the semiclassical prediction. We show in 
Fig-Elthe evolution of the correlations (Ax 2 ), and (Ap 2 ) 
of a trapped BEC with initial Uq = 5 when the nonlin- 
earity is tuned at time t = to U new — 0, 2, 7, and 
—0.5. It is noted that for U new = shown in Fig. [2a) 
the evolution of both position and momentum variances 
vary sinusoidally with equal amplitudes, consistent with 
the previous observation on the rotation of the error con- 
tour, while for U new ^ 0, the oscillation amplitudes are 
uneven, showing the effect of squeezing. 

This result clearly indicates that the amount of quan- 
tum squeezing can be controlled by controlling the mag- 
nitude of U new relative to Uq. From further simulations 
with different values of U new , it was found that, in gen- 
eral, with a decreasing U new < Uq such as Fig. [DJb) the 
minimum value that the position variance takes during 
the oscillatory dynamics decreases, indicating position 
squeezing while, with an increasing U new > Uq such as 
Fig- El^c), the maximum value that the position variance 
takes increases indicating momentum squeezing. This 
shows that the general effect of positive nonlinearity is 
to induce momentum squeezing. The case with a neg- 
ative Uq such as Fig. HJd) showed that the momentum 
variance quickly diverges as Uq becomes more negative. 

All these behaviors may be understood from the struc- 
ture of the extended Hamiltonian H ext and the changes in 
the infinite barrier potential term in Eq. (|16[) . It should 
be noted that the semiclassical predictions differ only in 
the amplitude for the U new = case while for U new ^ 0, 
there is also a phase difference so that over a long time, 
the semiclassical prediction gets worse. This feature due 



to the nonlinear effect agrees with the result observed in 
a different context of driven nonlinear systems [l^, and 
indicates that the semiclassical dynamics are better for 
predicting short time evolutions. 

Finally we consider the evolution of variances in an op- 
tical lattice instead of a harmonic trap. A salient aspect 
from the semiclassical equations is that, (x) and p are 
coupled in the presence of an optical lattice, unlike the 
case of a simple harmonic trap (a = (3 = 0) for which 
(x) and p are completely decoupled. Also, due to the 
e -p * factor, the effect of the optical lattice is signifi- 
cantly diminished when the wave packet is extended over 
the lattice i.e. when p 2 k 2 ^> 1 or p ^> Xol where Xol 
is the wavelength of the optical lattice. Likewise, when 
p is small, the effect of the optical lattice becomes more 
pronounced. 

Figure [3] shows what happens when one places a con- 
densate corresponding to Uq = 5 and position variance 
p 2 in an optical potential of varying wavelengths. It was 
found that for wavelengths Xol < p it is as if the optical 
lattice is not present, and shows ballistic expansion- like 
behavior. In such cases, the GPE simulation agreed very 
well with the semiclassical prediction. However, as Fig. 
SJa) shows, although the semiclassical equation captures 
the mean and variance it does not capture the detailed 
structure formed in the wave packet. For Xol — 2/5, we 
see in Fig. \%lb) that the semiclassical prediction starts 
to deviate from the calculation of the full GPE solution; 
this is understandable from Fig. 01b) where it shows 
the breakup of the wave function such that the Gaus- 
sian approximation cannot be valid (In fact, the Gaus- 
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FIG. 4: (Color online) Spatio-temporal plot of the proba- 
bility density |^(a;,i)| 2 obtained from the GPE when a BEC 
with the initial nonlinearity Uo — 5 and position variance p 2 
is placed in an optical lattice with wavelength (a) Xol = p, 
(b) Xol — 2p, (c) Xol — Wp, and (d) Xol = 20p. This shows 
that especially for Xol = p and 2p the Gaussian approxima- 
tion is not valid. 



sian approximation is not valid for the Xol = P case as 
well). For Xol = 10p and Xol — 20p shown in Fig.[3fc) 
and (d) the wavelength is sufficiently large that the wave 
packet gets "trapped" in one of the wells and since the 
initial state is not an eigenstate of the optical lattice the 
evolution in time was again oscillatory similar to that 
shown in Fig. O We found in the case of Xol ~> P 
that the smaller the value of the amplitude of the optical 
lattice a the greater is the position variance, indicating 
momentum squeezing; conversely, with a larger a one 
has increased momentum variance and position squeez- 
ing. A possible related phenomenon is the Superfluid- 
Mott insulator transition |21| (or a Josephson Junction as 
a two-site "lattice" , with the Rabi, Josephson, and Fock 
regimes in the order of increasing ratio of nonlinearity 
to tunnelling rate[I2].) To verify whether the Josephson 
Junction physics can be captured using semiclassical dy- 
namics, a new set of equations that assumes a Gaussian 



superposition state needs to be derived, and will be dealt 
with in a future work. 

The example with the optical lattice shows that there 
is rich physics to be explored, and that one needs to 
exercise caution in applying the semiclassical analysis 
to make detailed experimental predictions. The results 
show clearly that the quantitative differences between the 
semiclassical and the quantum result originate from the 
fact that the semiclassical limit is more or less a "coarse 
grained" representation of the full quantum dynamics. 
This is unavoidable since, fundamentally, classical me- 
chanics deal with point particles in phase space whereas 
quantum mechanics deal with wave functions. Indeed, 
the semiclassical approach such as the one discussed in 
this paper has been originally developed as one of the 
many attempts to bridge the gulf between the classical 
and quantum physics. 

In summary, we have derived a set of equations for 
atomic BECs which describe the condensate dynamics 
in the semiclassical limit. The semiclassical equations 
have been derived by applying the Gaussian squeezed 
coherent state ansatz for the condensate wave function. 
From the comparison of the semiclassical dynamics with 
the expected full quantum dynamics, it was found that 
the interatomic interactions tend to push the condensate 
deeper into the quantum regime. The advantage of these 
equations was that it was easy to get a more intuitive 
picture of the behavior of the system under considera- 
tion. Intuitive understanding could be obtained either 
from the classical Hamiltonian energy landscape argu- 
ment or from simply examining the equations carefully 
and considering the variables in various limits. Such 
an analysis is not possible with the GPE. Despite the 
additional advantages such as the very high numerical 
efficiency, if one wants to use the semiclassical equations 
as a predictive guide to the full quantum dynamics, one 
has to make sure that the Gaussian approximation is 
likely to be valid throughout the evolution. With this 
caveat, the semiclassical equations provide a quick diag- 
nostics and a verifiable "baseline" with which one can 
gauge the classical or quantum nature of a condensate 
dynamics under different situations. 

SC wishes to thank R. Reddy for discussions. 
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